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We study entanglement dynamics in a system consisting of a qubit dispersively coupled to a finite- 
temperature, dissipative, driven oscillator. We show that there are two generic ways to generate 
entanglement: one can entangle the qubit either with the phase or the amplitude of the oscillator's 
motion. Using an exact solution of the relevant quantum master equation, we study the robustness 
of both these kinds of entanglement against the effects of dissipation and temperature; in the limit 
of zero temperature (but finite damping), a simple analytic expression is derived for the logarithmic 
negativity. We also discuss how the generated entanglement may be detected via dephasing revivals, 
being mindful that revivals can occur even in the absence of any useful entanglement. Our results 
have relevance to quantum electromechanics, as well as to circuit QED systems. 

PACS numbers: 



There exists a long-standing interest in attempting to 
prepare and detect quantum states of macroscopic ob- 
jects or collective degrees of freedom. Such an experi- 
ment would be more than a mere "proof-of-principle" : 
it would provide a non-trivial test of our understanding 
of the quantum dissipative processes which cause such 
states to degrade with time, and thus enforce the quan- 
tum to classical transition. Recent advances suggest that 
submicron-scalc mechanical resonators could be excellent 
candidate systems in which to pursue this goal pp. Such 
resonators contain a truly macroscopic number of atoms; 
at the same time, they can be fabricated to have both 
high mechanical frequencies and quality factors. This 
suggests that one has some hope of cooling these systems 
to close to their quantum ground state, and that deco- 
herence due to the dissipative environment of resonator 
should be slow- in standard models, decoherence rates 
scale with the oscillator damping rate [2J. Nanomechan- 
ical resonators also have the advantage that they can be 
strongly coupled to (possibly coherent) electronic degrees 
of freedom; this has recently been demonstrated to allow 
sensitive position detection, approaching the fundamen- 
tal limits set by quantum back-action [3] . 

In this Letter, we analyze a promising electromechan- 
ical system where a dissipative mechanical resonator is 
dispersively coupled to a superconducting qubit: the 
state of the qubit simply shifts the frequency of the res- 
onator [HH]. Such a setup has the key advantage of being 
able to work with qubit states which are first-order in- 
sensitive to dephasing due to ever-present charge fluctua- 
tions [6j[7]. As such, the system is substantially different 
from the one analyzed in the seminal proposal of Armour 
et al. [5] , which made use of quickly-dephasing superposi- 
tions of charge states in the qubit. While the dispersive 
coupling allows for longer qubit coherence times, there 
is a price to pay: unlike the proposal of Ref. [HJ the two 
energy eigenstates of the qubit do not yield different aver- 
age forces on the oscillator. As such, generating entangle- 
ment is a slightly more involved affair. We demonstrate 



that there are two generic ways to generate non-classical, 
entangled oscillator-qubit states: one can entangle the 
qubit either with the amplitude of the oscillator's motion, 
or with its phase. We also study the robustness of these 
two kinds of entanglement against decoherence due to the 
dissipative environment, and discuss their detection us- 
ing coherence revivals. A fully analytical expression for 
the entanglement (as measured by the logarithmic nega- 
tivity) is obtained for the zero temperature case. In the 
finite temperature case, we make use of an exact solu- 
tion of the master equation to efficiently calculate the 
time-dependent entanglement. We stress that the sys- 
tem studied here has already been realized, both with 
nanomechanical resonators [S], and with superconduct- 
ing stripline resonators in circuit QED experiments |7J. 
Our study also sheds light on general questions of entan- 
glement dynamics in the presence of dissipation, driving 
and thermal noise. While entanglement has been stud- 
ied in a variety of qubit-plus-oscillator models [TUHH], we 
are not aware of a study involving the dispersive coupling 
considered here, nor of studies examining the effects of 
both driving and dissipation. 

System - We consider a damped mechanical oscillator 
(frequency um) which is dispersively coupled to a qubit 
(splitting frequency w q t>). Setting h = 1, the Hamiltonian 
takes the form H = H + H 1 with: 

H = (lj m + Aa z )(a 1 ! a + §) + + f(t) (a + a f ) 

(1) 

Here, A is the strength of the dispersive coupling, f(t) 
is an external force applied to the resonator, and H 1 de- 
scribes the damping of the oscillator (damping rate 7) 
by an equilibrium Ohmic bath at temperature T. We 
stress that such a dispersive coupling can be easily real- 
ized in systems having nanomechanical resonators cou- 
pled to superconducting qubits, as it emerges naturally 
from a Jaynes-Cumming type coupling in the relevant 
limit where w q b S> wm (i . Dispersive couplings have 
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also been achieved in the same way in recent circuit QED 
experiments coupling superconducting qubits to stripline 
resonators [12]. 

The dispersive coupling implies that the two qubit en- 
ergy eigenstates | |), | 1) each lead to different oscillator 
frequencies; equivalently, the effective frequency of the 
qubit depends on the energy of the oscillator. Unlike 
previous proposals (8], the two qubit states do not yield 
different oscillator forces, making entanglement genera- 
tion somewhat more subtle. One can easily show that if 
the oscillator starts in a thermal state and is not driven, 
then there is never any qubit-oscillator entanglement: the 
dispersive coupling simply leads to a statistical uncer- 
tainty in the qubit's frequency. Thus, to generate non- 
classical states, one must have a non-equilibrium state 
where (a(t)} is non-zero. To describe this, as well as the 
effects of dissipation and thermal noise, we will focus on 
the experimentally relevant regime where A, 7 <C u>m, 
and use the high-Q form of the Brownian-motion master 
equation for our system [J: 



P = ~i[Ho,p] +7(«oq- 



l)V[a]p + 7 n eq T>[tf}p 



(2) 



where for any operator A we define T>[A]p = ApA' — 
pA^A) /2. Here, 7 is the damping rate of the 
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resonator due the bath, Bose-Einstcin factor eval- 

uated at the bath temperature T and energy u>m, and 
r v is the intrinsic dephasing rate of the qubit. We con- 
sider the usual situation where the qubit energy relax- 
ation time is much longer than the dephasing time, and 
ignore T\ processes. 

As shown in Ref. [I] , Eq. ([2| can be solved for arbitrary 
T by using the fact that an initial Gaussian state remains 
Gaussian. With some work, this solution may be written 
in a physically transparent manner that is especially con- 
venient for entanglement calculations. For compactness, 
we focus on the relevant initial condition where at t = 
there is no qubit-oscillator entanglement, and the qubit 
state is \ip) = (| |) + | [))/^/2. Defining the displacement 
operator D\a] — exp[aa^ — a* a) and p atT ' as the partial 
trace Tr q b [/5|ct')((t|], the solution takes the form: 
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(l/2)D[ai(t)]-p eq [T\.&[ai(t)} 
(l/2)D[a l (t)].p eq [T\.&[a l (t)] 
[p iT (<)] f = (l/2)e^*Y(t) x 

D[S t (t)] • (p cq [T*{t)]e-^ t){fl+1 2 



(3a) 
(3b) 
(3c) 



Here, p eq (T) is the thermal equilibrium oscillator den- 
sity matrix at temperature T. As expected, py and 
PH are simply displaced thermal oscillator states, with 
the displacements a^it) denoting the conditional means 

^a(£)-|er)(o-y One has ay\ = —i (ujm ± A — 17/2) oty^ — 

if(t). The dynamics of py is more complex. We see 



that it also resembles a displaced thermal state. How- 
ever, the effective temperature T* is time dependent, 
and the dependence of the effective qubit frequency on 
the oscillator energy leads to a phase <fi(t). Defining 
a(t) = coth [(uj m /2T*) + i(f>/2], one finds: 



- 1 {a-{2n cq + l))-i\{a 2 -I) 



(4) 



with the initial condition T*(0) = T,0(O) = 0. In ad- 
dition, the displacements fi|^(t) arc generally not the 
expected, classical displacements ay^(t) appearing in 
Art/U' they only coincide at T = 0. One has: 

1 \ 1 (1 — Recr =F Hmo~)(a+ — a_) 
ay 1 = («±) ± T^TZ ( 5 ) 



2Re cr 



where we have defined 



-i (uj m ± ARe a) - - | u 



*f(t) (6) 



and 7 = 7 + AojA/Im a. For T = 0, a — > 1, and 
6ty^(t) = a + i_(t) = ay^(t). Finally, the parameter 
Y(t) describes the fact that due to the dissipative bath, 
the qubit-oscillator system will generally not be in a pure 
state ( Pn (T) = (n\p cq (T)\n)): 



Y(t) 



e Tipt exp 



iX f Q dt' (a + 2a + c 



(7) 



Entanglement- Similar to many recent works, we use 
the logarithmic negativity En to quantify the amount of 
qubit-oscillator entanglement in our system [T3J [H] . Em 
is a rigorous entanglement monotone applicable to mixed 
state systems, and is strictly zero for unentangled sys- 
tems. One has E^ — log 2 (2A/"+ 1), where the negativity 
M is the absolute value of the sum of the negative eigen- 
values of the partially-transposed density matrix p PT . 

To create non-classical, entangled qubit-oscillator 
states, we focus on the initial condition described above 
Eqs. Q: at t = there is no qubit-oscillator entangle- 
ment, and the qubit has been prepared (e.g. via a n/2 
pulse) in an equal superposition of its two eigenstates. 
The two qubit states in this superposition each lead (via 
the dispersive coupling) to different oscillator frequen- 
cies; this difference will be exploited to yield entangled 
states which, in the absence of dissipation or thermal 
noise, would have the simple form 



foMt))« [|t)KW) + ^ bt ll>KW>] 



(8) 



There are thus two generic ways to generate entangle- 
ment. The first is to entangle the qubit with the am- 
plitude of the oscillator's motion, i.e \(a)\. The second 
generic approach is to entangle the qubit with the phase 
of the resonator, i.e. arg(a). 

In the zero-temperaure limit (but still with 7 > 0) , the 
solution of the master equation in Eqs. (|3| yields an exact 
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expression for E^(t). Letting cos [0(t)] = \(a^(t)\ai(t))\ 
one finds: 



Af(t) 



(1 - |y|) 2 + 4|F|sin 2 6l- (1 - \Y\) 



(9) 



where for T = 0, Eq. for Y yields: 
\Y(t)\ 



e T ' pt exp 



2A/ dt' (|a T a;|sin</> n ) 
cos 6»(t) 



(10) 



with (f> n (t) 



arg 



As expected, Af (and 



a t (*)«!(*) 

hence the entanglement) is an increasing function of 
the distinguishability sin 2 6 of the two oscillator states. 
The factor \Y{t)\ describes as before the impurity of the 
oscillator-qubit system at times t > 0: as time progresses, 
the bath can distinguish the two oscillator states \a a ), 
thus reducing the oscillator-qubit entanglement. One 
might expect that \Y(t)\ should only depend on the his- 
tory of the overlap cos 6(t'); this is not the case. Instead, 
the decay of \Y(t) \ is sensitive to the history of sin </>-[■ j(f), 
i.e. the sine of the relative phase between the two oscil- 
lator coherent state amplitudes. Somewhat surprisingly, 
times when </>fj = it do not contribute to the decay of 
|F(t)|, even though the overlap between the two oscilla- 
tor states at such times is maximally small. 

In the case where T > 0, we again use our exact 
Eqs. ([3]) - (JsJ) to solve for the system dynamics. The form 
of p(t) given in Eqs. ^ then allows for a simple numeri- 
cal evaluation of Af: one simply converts p into a matrix 
in a basis of of displaced Fock states, and then numer- 
ically finds the partial transpose and the corresponding 
negative eigenvalues. For T > 0, E^{t) is not simply a 
function of the overlap (a^\ai)(t) and the purity Y(t) (as 
given by the full expression Eq. ([7]))- the entire matrix 
structure of p^i is relevant. 

Amplitude entanglement- To entangle the qubit with 
the amplitude of the oscillator's motion, we start at t = 
with the qubit in a superposition of its eigenstates, and 
the oscillator in a thermal state with (a) = 0. The oscilla- 
tor is then driven with a force f(t) — jctf cos[(wm + A)t] 
(af > 0). In the relevant limit of a high-Q oscillator 
where 7 -C A, f(t) will cause |a-|-(f)| to grow to a large 
value af, while |aj will be smaller by a large factor 
wm/t: hence, the amplitude of the oscillator's motion 
will become entangled with state of the qubit, and we 
would expect the qubit-oscillator entanglement to grow 
with time. However, at long enough times, the dissipa- 
tive bath coupled to the oscillator will destroy this en- 
tanglement: the bath can distinguish the two states \a a ), 
as described by Y(t). These two competing tendencies 
lead to being a non-monotonic function of time; this 
is shown in Fig. [l] Note that at T = 0, one can use 
Eqs. (J3j> and ^ to derive simple, exact expressions for 
Efj(t) for the amplitude entanglement setup; these will 
be presented elsewhere. 



1 

0.8 
0.6 1 
0.4 
0.2 




t (Y~ 1 ) 



FIG. 1: Log negativity En versus time. Blue lines correspond 
to the amplitude-entanglement setup, red lines correspond to 
the envelope of the oscillating En in the phase-entanglement 
setup. Solid lines are for n cq = 0, dot-dashed for n cq = 0.5, 
dashed for n cq = 1.0. For the phase-entanglement curves, we 
took ao = 0.76 to maximize the total integral of -Ejv(i) at 
n cq — 0. We also chose a/ = 3.74 so that the integral of 
-Ejv(t) for amplitude entanglement at n eq = is the same as 
the phase case. The inset shows the total time integral of En 
for both phase and amplitude entanglement as a function of 
n cq . In all cases, A = 0.01wm,7 = 10~ wm,^ = 0. 




FIG. 2: Left: Distance between effective oscillator displace- 
ments dtt(t),&i(t) entering the expression for pn. The ap- 
parent width of some curves is due to rapid oscillations with 
frequency \/ir. Right: Bose-Einstefn factor associated with 
the effective temperature T*(t) appearing in the expression 
for /5f j. Both cases are for the amplitude entanglement setup 
with af — 3.74. Temperatures are n cq = (black), n cq = 0.1 
(blue),n cq — 0.5 (green), n cq — 0.1 (red). A, 7 and Y v as in 



Fig. 1. For T > 0, one finds a CT 7^ a a 
the distinguishability of the a a . 



temperature reduces 



Fig. [T] also demonstrates that amplitude entanglement 
is dramatically suppressed by finite temperature. On a 
heuristic level, even though the two amplitudes a j , a j 
continue to have very different magnitudes at T > 
(they are of course independent of T) , the displacements 
<5|, ai which determine (c.f. Eq. (3c l)bccome less dis- 



tinguishable as T is increased. This behaviour is shown 
in Fig. [2] 

Phase Entanglement- To entangle the qubit with the 
resonator phase, we prepare the system at t = so that 
the qubit is again in a superposition of its two eigenstates, 
and the oscillator is in a state of motion characterized by 
the coherent state amplitude ao (e.g. one could drive the 
oscillator at lum, keeping the coupling to the qubit off 
until t = 0). At t = 0, we turn off all driving forces 
on the oscillator and let the coupled system evolve. The 
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FIG. 3: Upper panel: time-dependent log negativity En for 
the phase entanglement setup. Dashed line is for n oq = 0, 
dashed-dotted line for n cq = 1.0, solid for n cq — 5.0. Bottom 
panel: Qubit recoherences associated with phase entangle- 
ment, same parameters, qq = 0.76; X,y,T v as in Fig. 1. 



0.2 





n =0.1 

eq 




n =0.5 - 
eq 




n =1 . 

eq 




n =2 

eq 







e 1 1 1 1 1 

2 4 6 8 10 



a 2 Jn 
eq 

FIG. 4: Total entanglement (as measured by the time integral 
of En) for the phase entanglement setup, as a function of the 
initial coherent state amplitude q?o and n eq . For a finite n eq , 
there is an optimal value of qq which maximizes the total 
entanglement. In all cases, A = 0.01ojm,7 = 10~ ojm,^^ = 0. 



magnitudes of both coherent states a j , a j will be iden- 
tical, and will decay at a rate 7. In contrast, the phase 
of the oscillator coherent state will wind at a frequency 
determined by the qubit. We have thus prepared a state 
where the phase, not the amplitude, of the oscillator's 
motion is entangled with the qubit. 

At zero temperature, one can again obtain simple an- 
alytic expressions from Eq. ([9]); these will be presented 
elsewhere. Shown in the top panel of Figj3]is E^(t) for 
the phase entanglement setup. The entanglement drops 
to zero periodically at a frequency A/7r: at these times, 
the phases of the two oscillator states aj, aj are aligned, 
implying cos 6=1. This alignment also implies that 
the coherence of the qubit is restored, as is shown in the 
bottom panel of Fig. [3j where we plot 2|Tr/5-j-jJ. 

It is interesting to compare the phase and amplitude 
entanglement setups; this is done in Fig. [l] where we 
plot only the enveloped of the oscillating entanglement 
in the phase case. As can be seen from the inset, phase 
entanglement is more robust against non-zero T than 
amplitude entanglement; the total integral of Ej^(t) de- 
cays far more slowly as a function of n eq in the phase 
case, for parameters that yield the same entanglement 
at n cq = 0. Further, while both kinds of entanglement 
are suppressed by non-zero T, entanglement in the phase 
case can remain large for short times (i.e. for times t such 
that 1/A <C t <C 1/7). Thus, phase entanglement shows 
a certain increased resilience against thermal dissipation 
compared to amplitude entanglement. For a fixed bath 
temperature, there is in general an optimal value of ao 
which maximizes phase entanglement; this is shown in 
Fig. [4] For ao = 0, E]y(t) = for all times, where as for 
too large an ao, the bath very rapidly distinguishes the 
two oscillator states in the superposition. 

Detecting Entanglement - We now turn to the detec- 
tion of phase entanglement using coherence revivals. As 
shown in Fig. [3] such entanglement leads to qubit reco- 
herences: the magnitude of the qubit 's off-diagonal den- 
sity matrix element |Tr/Sjj(i)| is non-monotonic in time. 
This quantity represents the time-dependent dephasing 
of the qubit, and is measurable via either a standard 



Ramsey interference experiment [15 , or via state tomog- 
raphy. While such revivals of coherence have been used 
to detect non-classical states in other situations [16] , and 
have been proposed as a way to detect entanglement in 
NEMS [8 , one must be careful: it is possible to have 
coherence revivals without any qubit-oscillator entangle- 
ment. In our system, a purely thermal state oscillator 
state with (a(t)) = yields dephasing revivals, but zero 
qubit-oscillator entanglement. 

Despite this caveat, one can still use dephasing revivals 
as a proxy for detecting entanglement. The basic idea is 
that the Fourier spectrum of the time-dependent dephas- 
ing lets one unambiguously distinguish an initial thermal 
state with a = and E^it) = 0, from a phase-entangled 
state where a^0. This difference is essentially the num- 
ber splitting effect discussed in Refs. |3J HZJ EH] 7 and mea- 
sured in Ref. [T^]: the peaks in the Fourier spectra are 
directly related to the number distribution in the oscil- 
lator. Thus, in a thermal state, one expects the peaks 
to follow a simple geometric series. In contrast, when ao 
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FIG. 5: S(u), the real part of the Fourier transform of the 
qubit 's time-dependent coherence 2 |Tr/5f|(i)|, for the phase- 
entanglement setup. The red-dashed curve is for n cq = 0.5, 
ao — 0. There are revivals of coherence, but strictly zero 
qubit-oscillator entanglement; S(u)) show peaks whose area 
follows a geometric series. In contrast, the solid blue curve 
corresponds to ao = 1.23, n cq = 0.5. Here, one gets both co- 
herence revivals and entanglement. The dephasing spectrum 
is markedly different: the peak areas do not form a geometric 
series. The inset shows the integrated area for each spectra in 
the two cases. A = O.OIwm, 7 = 10 _5 um and T v — 10 _3 cjm- 



is non-zero, the distribution begins to resemble more the 
Poisson distribution associated with a number state. Ob- 
serving this difference, and comparing it to theory, would 
be a convincing way to detect the phase entanglement we 
have described. The difference between the two spectra 
is shown in Fig. [5] 

Conclusion- We have studied entanglement in a dis- 
persive qubit - oscillator system, identifying two generic 
kinds of entanglement: phase and amplitude entangle- 
ment. We have shown that in general, phase entangle- 
ment is more robust against the effects of dissipation, and 
have shown how it may be detected through the analysis 
of revivals in the qubit's coherence. 
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